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We have studied theoretically the basic operation of a 
quantum feedback loop designed to maintain a desired phase 
of quantum coherent oscillations in a single solid-state qubit. 
The degree of oscillations synchronization with external har- 
monic signal is calculated as a function of feedback strength, 
taking into account available bandwidth and coupling to envi- 
ronment. The feedback can efficiently suppress the dephasing 
of oscillations if the qubit coupling to the detector is stronger 
than coupling to environment. 

The principle of feedback control is used in a wide va- 
riety of physical and engineering problems. In particular, 
it can be applied in a straightforward way to tune the os- 
cillation phase of a harmonic oscillator in order to achieve 
a desired synchronization with some reference oscillator. 
An intriguing and fundamental question is whether con- 
tinuous feedback can be used to control quantum sys- 
tems; for instance, if it is possible or not to tune the 
phase of quantum coherent (Rabi) oscillations in a qubit 
(two- level system). 

At first sight the quantum feedback seems to be im- 
possible because according to the "orthodox" collapse 
postulates the quantum state is abruptly destroyed by 
the act of measurement. However, jas was shown two 
decades ago, in particular by Leggett,El in a typical solid- 
state setup the collapse of a qubit state should be consid- 
ered as a continuous process rather than as instantaneous 
event. The reason is typically weak coupling between the 
quantum system and the detector and also the finite noise 
of the detector, so that it takes some time until accept- 
able signal-to-noise ratio is reached and the measurement 
can be regarded as completed. 

While the Leggett's theory as well as the majority 
of similar approaches can describe only ensembles of 
quantum systems, the theory describing the gradual col- 
lapse oi-^i-single solid-state qubit was developed only 
recentlyJjp (A similar problem in quantum optics was 
solved much earlier - see, e.g. Refs.l3^Ll and references 
inO.) Basically, the theory says that the evolution of a 
single quantum system due to continuous measurement is 
governed by the information continuously acquired from 
the detector. Similarly to classical probability, the Bayes 
formulala which naturally takes into account incomplete 
information from the detector, can still be applied to the 
density matrix of the measured |-Quantum system; thus 
the formalism is called Bayesian.u 

In case of a poor detector the extra noise acting back 
onto the input disturbs the measured system stronger 
than the limit determined by the uncertainty principle; 
this leads to gradual decoherence of the measured sys- 
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FIG. 1. Schematic of the quantum feedback loop main- 
taining the quantum oscillations in a qubit. 



tem. In contrast, when measured with a good (quantum- 
limited) detector, the quantum system does not loose the 
coherence (even though the quantum state evolves ran- 
domly)' moreover, its density matrix can be gradually 
purifiedd that basically means acquiring as much infor- 
mation about the system as permitted by quantum me- 
chanics. 

Since the Bayesian formalism allows us to monitor the 
continuous evolution of a quantum system in a process 
of measurement, this naturally gives rise to a possibility 
of continuous feedback control of a quantum system. In 
this paper we wilLstudy the operation of a feedback loop 
proposed in Ref.Q and designed to maintain a desired 
phase of quantum coherent oscillations in a solid-state 
qubit. (Quantum feedback in quantum opticajijas-been 
proposed and studied earlier - see, e.g., Refs.Qclii3.) In 
particular, we will study the dependence of the loop op- 
eration on the feedback strength, available bandwidth, 
and dephasing due to environment. 

As an example of the measurement setup (Fig. |l|) we 
consider a qubit represented by a single electron in a dou- 
ble quantum dot (DQD), the location of which is mea- 
sured by a qiiantum point contact (QPC) nearby in a way 
used in Refill. If the electron is in the dot 2 (state |2)) 
which is closer to QPC than dot 1, then the QPC tunnel 
barrier is higher and so the average current I2 through 
QPC is smaller than the average current Ii correspond- 
ing to the electron in the dot 1 (state |1)). Consequently, 
from the QPC current one gets information about the 
electron location. We consider a realistic case of weak 
response, A/ = Ji - /2 < /o = (/i + h)/'^.- In this 
case the measurement time 6*7/2 (A/)^, which is neces- 
sary to achieve signal-to-noise ratio equal to 1 (here Si 
is the shot noise of the QPC current), is much larger 
than e//o, so the QPC current I{t) is continuous on the 
measurement timescale and we do not need to consider 
individual tunneling events in QPC. 

The evolution of the qubit density matrix p during the 
measurement process ia pijiescribcd within the Bayesian 
formalism by equationsl 
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Pii = -P22 = -2 Impi2 + P11P22 -5— [I[t) - -fo], (1) 
n oj 

. £ . H , 

P12 = 1 ^ P12 + 1 T (Pii ~ P22 J 
n n 

- (pii - P22)^ [/(i) - Iq] P12 - 7P12 , (2) 

where e and H are, respectively, the energy asymmetry 
and tunneling strength of the qubit [the qubit Hamil- 
tonian is Hqb = {e/2){c\c2 - c\ci) + H{c\c2 + c\ci)\, 
and 7 = 7d + 7e is the dephasing rate due to the detec- 
tor uonideality (7^;) and coupling with the environment 
(7e).EZl Theoretically, 7,; = when qubit is measured 
by a QPC; however, if instead of QPC we use a single- 
electron transistep-|(SET), then dephasing is usually 
quite significantQila (except tliejCase when the SET op- 
erates in a cotunneling regimellj'EHl) . 

Notice that the ensemble dephasing rate F = 7 -|- 
(A/)^/4S'/ is larger than 7 because of different evolution 
of the ensemble members due to random I{t). Ipjdividual 
realizations can be simulated using the formulaa 



/(t)-/o-(pii-p22)A//2 + eW, 



(3) 



where (^{t) is the pure white noise with spectral density 
5^ = Si. If Eqs. (0)-(§) are averaged over ^(i) (we use 
Stratonovich definition for stochastic differential equa- 
tions), then we get usual ensemble-averaged equations 
for qubit evolution (terms proportional to A/ will disap- 
pear and 7 will be replaced by F). 

It is natural to characterize the effect of extra de- 
phasing 7d by the detector idealitjii| (efficiency) rj = 
I/[l-h7d45'7/(A/)2]. One can showfclthat r] = lh/2edf 
where is the total energy sensitivity of the detector 
[td = (eiEo)^''^ where Co is the usual (output) energy sen- 
sitivity and ei is a similar quantity characterizing back- 
action to the input]. So, an ideal case r/ — I corresponds 
to a detector with quantum-limited sensitivity. 

To realize a feedback loop (Fig. we can moni- 
tor the qubit evolution using the detector current I{t) 
plugged into Eqs. (|l|)-(||). Then the qubit state is com- 
pared with the desired state, and the difference sig- 
nal is used to control the qubit parameters H and/or 
e. In the example studied in this paper the feedback 
loop is designed to stabilize the quantum oscillations 
of the state of a symmetric qubit (e = 0), so the de- 
sired evolution is pii{t) = 1 — P22{t) = [1 + cos(f2t)]/2, 
Pi2{t) = P2i{t) = isin(rit)/2, where the frequency is 
n = (4^2 + £2ji/2/7j ^ 2H/h. As a difference ("error") 
signal we use the phase difference Aip (| A0| < tt) between 
the desired value 0o (t) — (mod 27r) and the moni- 
tored value (j){t) = arctan{2Impi2(t)/[jOii(t) — P22{t)]}- 
This difference is used to control the qubit parameter H 
(changing the barrier height of DQD); here we study a 
linear control: Hft, = (1 — Fj-x A(j))H, where F is the 
dimensionless feedback factor .E3 

In this paper we neglect additional time delaycl in the 
feedback network, however, we take into account the fi- 
nite bandwidth of a line carrying detector current (that 
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FIG. 2. Correlation function K^ij) of the qubit quantum 
oscillations for C = 1 and feedback factors F = Q (thin solid 
line), 0.05 (thick solid line), and 0.5 (dashed line). Nondecay- 
ing oscillations are due to synchronization by the feedback. 



is a critical parameter for a possible experiment). More 
specifically, we average the current I{t) with a rectan- 
gular window of duration Ta, Ia{t) = J^_^ I{t')dt', 
before plugging it into Eqs. (jl|)-(^), so that the "avail- 
able" density matrix Pa{t) differs from the "true" density 
matrix p{t). Also, to compensate for the corresponding 
implicit time delay, we use A0 = (j)a — il,(t — kTo) with 
K — 1/2 (we tried various n and found that k = 1/2 
provides the best operation of the feedback loop) . 

Let us start with the case of ideal detector, -q — 1, 
absence of extra environment, 7e = 0, and infinite 
bandwidth, Tq = 0. Figure |^ shows numerically calcu- 
lated correlation function Kz{t) = {z{t + T)z{t)) where 
z = pii — P22, for several feedback factors: F = 0, 
0.05, and 0.5. Jie curves are obtained using Monte 
Carlo simulatiorou of the measurement process for mod- 
erately weak coupling between the qubit and detector: 
C = n{AD3/SiH = I (notice that the Q-factor of 
oscillationscJ is equal to 8/C, so C = I is still a weak 
coupHng). In absence of feedback {F — 0) the corre- 
lation function decays to zero (Fig. |^) while for finite 
feedback factor the correlations remain for indefinitely 
long time (of course, assuming perfect reference oscilla- 
tor which determines the desired evolution). The non- 
decaying correlations show that the quantum feedback 
loop really provides the synchronization of quantum os- 
cillations. The degree of synchronization depends on the 
feedback factor F. One can see that for a moderate value 
of -F = 0.5 the synchronization is already very good [the 
ideal case would be Kz{t) = cos(fiT)/2]. 

For analytical analysis we take into account-that in the 
ideal case 7d = 7e = the qubit state is pure,Q and using 
Eqs. (0)-(p|) start with the equation 
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A(j) ■ 



AI AI ^ \ 2FH ^ ^ 
sm — I — cos -I- M ^ A(j), (4) 
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which assumes the absence of 2tt phase slips (good or 
moderate synchronization). For weak coupling (C/8 <C 
I) we can neglect the first term in parentheses and av- 
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FIG. 3. Dots: asymptotic amplitude Akz of Kz{t) oscilla- 
tions as a function of feedback factor F for several couplings 
with the detector, C = 0.5, 1, and 2. Solid lines: analytical 
approximation Akz = exp(--C/16i^)/2. Dashed lines: corre- 
sponding numerical results for /2. 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
F (feedback factor) 

FIG. 4. Synchronization degree D as a function of feedback 
factor F for several values Ta of detector signal averaging: 
Ta/T = 0, 1/3, and 2/3, where T = 2n/D,. Dashed line 
D = exp(— C/32f) almost coincides with the upper curve. 
Dotted line corresponds to "direct" feedback with Ta = T/IO. 



erage the random term over sin^ assuming almost har- 
monic evolution that leads to the simplified equation 



d ^ ^ ~ 2FH ^ ^ 



(5) 



where ^(t) is the white noise with spectral density — 
{AI)'^/2S. This equation describes a particle diffusion 
in the parabolic potential (we again assume |A0| < tt). 
The corresponding Fokker-Planck equation has an exact 
solution which is used to calculate the correlation func- 
tion K^{t) « (cos[A0(i) - A(j}{t + t)]) cos nT/2. In this 
way we obtain the analytical expression 
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which fits well the Monte-Carlo results when C/8 ^ 1 
and C/IQF < 1 (weak coupling and moderate or good 
synchronization). As an example, the dots in Fig. ^show 
the numerically calculated (using the least-mean-square 
fit) asymptotic amplitude Akz of Kz{t) oscillations (at 
T — > oo) as a function of the feedback factor F for three 
values of the coupling C, while solid lines show the cor- 
responding analytical curves Akz = exp(— C/16F)/2. 

The correlation function Ki{t) = {I{t + t)/(<)) of the 
detector current I{t) is somewhat similar to-Kz{T), how- 
ever, it also has the decaying contributions due to cor- 
relation Kz$^ and a (5-function contribution due to the 
detector noise. The analytical result for the same regime 
as above, 
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also agrees well with the Monte Carlo results. 



The spectral density Si{u!) of the detector current can 
be obtained as a Fourier transform of Ki{t). While in 
absence of feedback the quantum oscillations in the qubit 
can provide only a moderate peak of 5/(w) around fre- 
quency f2 (the peak height cannot be larger than 4 times 
the noise pedestaled) , the feedback synchronization leads 
to the appearance of a (5-function at the frequency of de- 
sired oscillations. (In principle the desired frequency can 
differ a little from $1; however, in this case the perfor- 
mance of the feedback loop worsens.) 

Besides the correlation function and spectral density, 
we have studied one more characteristic, D, of the syn- 
chronization degree. We define D as the average scalar 
product of the unity-length vector on the Bloch sphere 
corresponding to the desired state and the vector corre- 
sponding to the actual state of the qubit. The equivalent 
definition is D = 2{TTppii) — 1, where pd is the density 
matrix of the desired pure state. [The so-called fidelity 
is equal to eithei-^D -I- l)/2 or {D + l)/2, depending 
on the definitionJlj] Perfect synchronization corresponds 
to D = 1. It is simple to show that in the limit of weak 
coupling and for symmetric distribution of Ac/) (unshifted 
desired frequency), Akz coincides with D^/2. Notice, 
however, that at moderate coupling, D'^/2 (see dashed 
lines in Fig. ^) is significantly closer to the analytical 
result than Akz- 

Upper solid line in Fig. ^ shows the dependence of D 
on the feedback factor F for C = 1 and Tq = 0. One 
can see that D is proportional to F for small F ("soft" 
onset of the synchronization) and D is asymptotically 
approaching 1 at large F. The analytical result D = 
exp(-C/32F) (dashed line in Fig. |) is very close to the 
numerical results at moderate and good synchronization. 

Finite available bandwidth of the detector current I{t) 
(finite averaging time Tq in our formalism) worsens the 
performance of the quantum feedback loop. The solid 
lines in Fig. show the dependence of the synchroniza- 
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FIG. 5. Dependence D{F) for C = 1, Ta = 0, and several 
magnitudes of dephasing due to environment; — 0, 0.1, 
and 0.5. Dashed and dotted lines correspond to de = and 
limitation of Hfb by and H/2, respectively. 

tion degree D{F) for Ta/T = 0, 1/3, and 2/3, where 
T = 27r/il is the oscillation period. Obviously, a sig- 
nificant information loss occurs when Ta becomes com- 
parable to T, leading to a decrease of D. The curves 
D{F) saturate at large F allowing us to introduce the 
dependence Dmax{T)- Calculations for the parameters of 
Fig. H show pretty good synchronization, Dmax = 0.993, 
for Ta = T/30, while Dmax = 0.98, 0.92, and 0.57 for 
Ta = T/10, r/3, and 2T/3, respectively. 

The main potential practical importance of the quan- 
tum feedback is the ability to suppress the effect of the 
qubit dephasing caused by interaction with the environ- 
ment (see Fig. |l|). This can be used, for example, for 
qubit initialization in a solid-state quantum computer. 
Solid lines in Fig. ^ show the dependence D{F) for sev- 
eral magnitudes of the dephasing due to environment, 
de = 0, 0.1, and 0.5, where 4 = 7e/[(A/)2/4S'/] is the 
ratio between the qubit coupling to the environment and 
to the detector (we still assume an ideal detector). First 
of all, we see that the feedback still maintains the qubit 
phase synchronization for infinitely long time. However, 
for finite de the degree of synchronization D saturates 
at a level less than unity. We have studied numerically 
the dependence Dmax{de) for C = 1/2, 1, and 2 (while 
Ta = and T] — \) and found a linear dependence at 
small de- D„iax — 1 — 0.5 d^. [A little better formula 
Dmax — 1 — 0.5 de/ (1 -l- de) works reasonably well up to 
de ^ 1.] This means that the feedback loop can effi- 
ciently suppress the qubit dephasing due to the coupling 
to the environment if this coupling is much weaker than 
the qubit coupling to a nearly ideal detector. 

Notice that the solid lines shown in Figs. ^ and |^ are 
calculated assuming the feedback control of the tunnel 
matrix element Hfi, — H[l — F x Ac/)] even when Hfi, 
becomes negative (this is also an assumption for the ana- 
lytical results). To eliminate this unphysical assumption 
we have also performed numerical calculations with re- 
strictions Hfb > and Hfi, > H/2. This leads to rather 



minor modifications of the presented curves (dashed and 
dotted lines in Fig. H show the results for de — and 
Ta = 0). However, important difference is that D{F) 
goes down at large F, so the optimum Dmax is achieved 
at some finite value of F. More detailed study of this 
problem will be presented elsewhere. 

Besides the discussed feedback based on Acj) calcula- 
tion, we have also studied a "direct" feedback loop in 
which Hfbit)/H - 1 = F{2[Ia{t) - Io]/AI - cos[17(i - 
Tq/2)]} sin[il{t — Ta/2)] (we call it also a "naive" feedback 
because this control formula is easily designed from the 
naive assumption that the detector current directly fol- 
lows the evolution of pii). Direct feedback is much sim- 
pler for experimental realization since it does not require 
real-time solution of the Bayesian equations (directi'eed- 
back in quantum optics has been studied in Refs.EHl3 EJ). 
Quite surprisingly for us, the direct feedback can also 
provide a good phase synchronization of quantum oscil- 
lations if F/C is close to 1/4 (see dotted line in Fig. U). 
However, it requires more careful choice of F and than 
for the Bayesian feedback, and also suffers more signifi- 
cantly from the restriction on Hft variation. The results 
in more detail will be discussed elsewhere. 

Experimentally, besides the realization of quantum 
feedback control of a DQD continuously measured by 
a QPC, one can also think about the qubit based on 
a single- Cooper-pair box measured by a single-electron 
transistor (see discussion iro). This realization can be 
preferable because of a rapid progress of metallic single- 
electronics technology. However, the problems are high 
output impedance of the single-electron transistor and its 
nonideality as a quantum detector. The third potential 
realization can be based on SQUIDs. For any realiza- 
tion the major problem is the bandwidth: the feedback 
should be at least faster than the qubit dephasing. Be- 
cause of that, the quantum feedback of a solid-state qubit 
should probably be attempted only after the realization 
of recently jproposed Bell-type two-detector correlation 
experiment JHj which would show the possibility of quan- 
tum monitoring, the first step to quantum feedback con- 
trol. 
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